Late time acceleration in Brans Dicke Cosmology 



S. Sen* and A. A. Sent 
Harish Chandra Research Institute, Chhatnag Road, Jhusi. Allahabad 211 019 India 

(February 7, 2008) 

In this work we have investigated the possibility of having a late time accelerated phase of the universe, sug- 
gested by recent supernova observation, in the context of Brans Dicke (BD) theory with potential having a time 
dependent mass squared term which has recently become negative and a matter field. We find that while a per- 
fect fluid (pressureless and with pressure) cannot support this acceleration, a fluid with dissipative pressure can 
drive this late time acceleration for a simple power law expansion of the universe. We have also calculated some 
cosmological parameters in our model to match with observations. 
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I. INTRODUCTION 

A lot of activity has been triggered by two recent obser- 
vations [pl^|] on the explosion of type la Supernovae. These 
datas favour the existence of a new kind of matter with posi- 
tive energy density dominant at present universe and is also 
responsible for the present acceleration of the universe ac- 
counted by its negative pressure. This along with the observed 
location of the first acoustic peak of CMB temperature fluctua- 
tion corroborated by the latest BOOMERANG and MAXIMA 
data [ppj], favours a spatially flat universe whose energy den- 
sity is dominated by a cosmological constant like term. Obvi- 
ously the first natural choice to represent such special matter 
was the cosmological constant A For a flat matter dom- 
inated universe with A having Ha ~ 0.72 in Einstein gravity 
best fits the data sets. But the candidature of A as the con- 
stituent of the major energy density is troubled by the fact that 
it has an energy scale which is ~ 10~ 123 lower than normal 
energy scale predicted by the most particle physics models. So 
to find some alternative candidate for this acceleration a dy- 
namical A |Q| in the form of a scalar field with some self inter- 
acting potential Q is considered whose slowly varying energy 
density mimics an effective cosmological constant. The idea 
of this candidate, called quintessence j7|, is borrowed from 
the inflationary phase of the early universe, with the differ- 
ence that it evolves at a much lower energy scale. The energy 
density of this field, though dominant at present epoch, must 
remain subdominant at very early stage and has to evolve in 
such a way that it becomes comparable with the matter den- 
sity n m now. This type of specific evolution, better known as 
"cosmic coincidence" ^ problem, needs several constraints 
and fine tuning of parameters for the potential used to model 
quintessence with minimally coupled scalar field. A new form 
of quintessence field called "tracker field" jiot ] has been pro- 
posed to solve the cosmic coincidence problem. It has an 
equation of motion with an attractor like solution in a sense 
that for a wide range of initial conditions the equation of mo- 
tion converges to the same solution. 

There are a number of quintessence models which have 
been put forward and most of which involve minimally cou- 



pled scalar field with different potentials dominating over the 
kinetic energy of the field. A purely exponential potential 
is one of the widely studied cases [11]. Inspite of the other 



advantages the energy density is not enough to make up for 
the missing part. Inverse power law is the other potential ( 
[|lO[]) that has been studied extensively for quintessence 
models, particularly for solving the cosmic coincidence prob- 
lem. Though the problems are resolved successfully with 
this potential, the predicted value for the equation of state 
for the quintessence field, jq, is not in good agreement 
with the observed results. In search of proper potentials that 
would eliminate theproblems, new types of potentials, like 
V [coshX(f) - 1]P and V sin h(ay/k [|@] have 
been considered, which have asymptotic forms like the inverse 
power law or exponential ones. Different physical consider- 
ations have lead to the study of other types of the potentials 
also [|l6|]. Recently Saini et al [|l7]] have reconstructed the po- 
tential in context of general relativity and minimally coupled 
quintessence field from the expression of the luminosity dis- 
tance dh (z) as function of redshift obtained from the observa- 
tional data. However, none of these potentials are entirely free 
of problems. Hence, there is still a need to identify appropri- 
ate potentials to explain current observations [pTj], 



Most of the studies regarding accelerated expansion have 
been done with a minimally coupled scalar field representing 
the quintessence. It has been recently shown by Pietro and 
Demaret Jl2] ] that for constant scalar field equation of state, 
which is a good approximation for a tracker field solutions, the 
field equations and the conservation equations strongly con- 
strain the scalar field potential, and most of the widely used 
potential for quintessence, such as inverse power law one, ex- 
ponential or the cosine form, are incompatible with these con- 
straints. The minimally coupled self interacting models will 
also be ruled out if the observations predict that the missing 
component of the energy density obeys an equation of state 
p = 7/9 with 7 < — l(p > 0) , and these sort of equation 
of state is in reasonable agreement with different observations 
[|l|]. Also the inequality dH 2 {z)/dz > ZH mQ H {l + z) 2 
should satisfy for minimally coupled scalar field and its vi- 
olation will certainly point towards a theory of non Einstein 
gravity such as scalar tensor theories where the scalar field is 
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non minimally coupled to gravity. 

There have been quite a few attempts of treating this prob- 
lem with the non-minimally coupled scalar fields. Scal- 
ing attractor solutions are available in the literature with the 
exponential [ jl6[ l and power law [16, 181] potentials in non- 
minimally coupled theories. Faraoni [|19|] have studied dif- 

j 2 

ferent potentials with a non-minimal coupling term ipR^k- for 
the present acceleration. There have been different approach 
also for solving the problem in general scalar tensor theory, 
sometimes called extended or generalised quintessence, not 
only because this theory is considered to be the most natural 
alternative to general relativity, there are other strong motiva- 
tions j20| ] also. People like Bertolo et al [|T|, Bertolami et al 
[p2[], Ritis et al [ p3[ ] have found tracking solutions in scalar 
tensor theories with different types of power law potential. In 
another work Sen et al [^4j have found the potential relevant 
to power law expansion in Brans Dicke cosmology. Like Saini 
et al jr7|], Boisseau et al [ p5[ ] have reconstructed the potential 
from the luminosity-redshift relation available from the obser- 
vations in context of scalar tensor theory. 



Very recently McDonald [|30|] has investigated the possibil- 
ity of modelling a dynamical cosmological constant with a 
scalar field which has undergone a very recent phase transi- 
tion. For this he has considered a standard </> 4 potential for the 
scalar field with an additional time dependent mass squared 
term in the potential which has become negative very recently. 
For this kind of model, phase transition occurs very recently 
atredshiftz < 1.2. 

In this paper we have investigated whether non minimally 
coupled self interacting scalar field such as a Brans-Dicke 
(BD) type scalar field with this kind of potential can success- 
fully drive the late time acceleration for the flat universe. In 
the context of Brans Dicke(BD) theory [^9|] with a self inter- 
acting potential and a matter field, the action is given by 



S = 



(1) 



(We have chosen the unit 8ttGq — c = 1.) 
We have chosen the cf) 4 potential with a time dependent mass 
squared term which has already become negative after a phase 
transition in recent time [BQ] : 



where 



M 2 (*)=Mo 2 



Mo. 
R n 



(2) 



(3) 



Mo = fio 2 Rc n - A is a constant and n is an integer. The time de- 
pendent mass squared term with integer n can arise naturally 
in plausible models and one can have a detailed discussion in 



As a matter field we would consider first perfect fluid and 
then a fluid having negative pressure. An effective negative 
pressure and hence an acceleration can be achieved by dis- 
sipative mechanism modelled commonly by fluid viscosities. 
It has been proposed recently that the CDM must self inter- 
act in order to explain the detailed structure of the galactic 
halos fl27[]. This self interaction will naturally create a vis- 
cous pressure whose magnitude will depend on the mean free 
path of the CDM particles. An effective negative pressure in 
CDM can also be created from cosmic anti friction which is 
closely related to particle production out of gravitational field 
fl28]]. Since the negative pressure can be modelled in different 
ways, we are not apriori assuming any specific model for this 
negative pressure. 

In this work we find that it is not possible to have a late 
time power law accelerated expansion when the CDM is a 
perfect fluid, but a dissipative CDM fluid in BD cosmology 
with such a potential like (^) can successfully drive a late time 
accelerated expansion. In the next section we treat the field 
equations and find the solutions in both the cases. We also 
calculate some cosmological parameters to match the acceler- 
ated model with observation. The third section is the conclud- 
ing section where we have discussed different features of this 
model. 



II. FIELD EQUATIONS AND SOLUTIONS 

The gravitational field equations derived from the action ([!]) 
by varying the action with respect to the metric is, 

+t[<^ ; „ - 9»» a <P] ~ 9»J^- (4) 

where T^ v represents the energy momentum tensor of the mat- 
ter field. We have assumed the matter content of the universe 
to be composed of a fluid represented by the energy momen- 
tum tensor 



T^u = {p + P)v fl v v + Pg^u, 



(5) 



where p and P are the energy density and effective pressure of 
the fluid respectively and u M is the four velocity of the fluid i.e, 
v^v^ = —1. The effective pressure of the fluid includes the 
thermodynamic pressure p and a negative pressure ir, which 
could arise either because of the viscous effect or due to par- 
ticle production, i.e, 



P = p + 7T 



(6) 



The wave equation that follows from equation (|l|), by varying 
the action with respect to the scalar field <j> is 



T 
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2ui + 3 2cu + 3 
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2V(4>) (7) 
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For our choice of potential (g), the field equations (Q) and the 
wave equation (0) becomes respectively 
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(8) 



(p + ir) 



(9) 



and 



i? • = p - 3(p + tt) _ 1 
ii 2u + 3 2cj + 3 



2A</> 4 



(10) 



We have assumed standard Friedman-Robertson- Walker met- 
ric with the signature convention (—,+,+,+) and R is the 
scale factor. We restrict ourselves for spatially fiat metric only. 
We work in Jordan frame. One interesting thing about BD the- 
ory in Jordan frame is that the conservation equation holds for 
matter and scalar field separately. Or in a slightly different 
way, the Bianchi Identity along with the wave equation (0) 
gives the matter conservation equation 



3|(, 



P 



tt)=0 



(ID 



We assume that both the scale factor and the scalar field evolve 
as the power function of time 



/»> = ./?„ [ ~ 
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and d> = 



[ r 



(12) 



where the subscript refers to the values of the parameters at 
the present epoch. In order to get accelerated expansion for 
such a evolution of the universe the deceleration parameter 
has to be negative, which immediately restricts the parameter 
a to be greater than 1 . For such an expansion the solution for 
the matter density is 



(3-2 



where 



3ac 



t 



Pet 



2a + (3(l + a)-(3 2 (l+cu) 



2-/? 



(13) 



(14) 



First we consider normal perfect fluid with no negative pres- 
sure i.e, 7r = in equation (||). Then, for the power law evo- 
lution, the thermodynamic pressure of the fluid becomes 



P = P, 



where 



Pc = 



(2-/3- 3a)d> 
t 

I O 



2a + /3(l + a)-/3 2 (l+cj) 



2-/3 



(15) 



(16) 



Power law solution is consistent with the field equations 
(g) and © only if 



and an — (5 = 2 i.e, an 



4 
3 



(17) 



So the acceleration demands n < |, and as n is an positive 
integer, it restricts n = 1 in our model. Again from equation 
(fijl) the weak energy condition (p > 0) demands 



u> < 3a 



(18) 



From equation (jlj) and ( |l3| ) it is clear the perfect fluid follows 
an equation of state of the form p = j m p- The index 
given by 



is 



7n 



2-/3 
3a 



1. 



(19) 



where j m lies within the interval < 7 m < 1. This restricts 
a within the range | < a < |. Infact for present matter 
dominated universe (j m = 0), a = |. But this does not 
satisfy the criteria for acceleration (a > 1) and hence for 
a simple power law type expansion the universe decelerates 
with a perfect fluid CDM (0 < j m < 1) with a potential (||) 
in BD theory. 

Now we consider a CDM which has a dissipative effect and 
we are particularly interested in a present day universe i.e, 
p = 0. Under such condition equation (p4[) takes the form 



3|(p + 7r)=0 



(20) 



As is mentioned earlier this type of dissipative effect in FRW 
cosmology can be modelled in two ways. Generally the dissi- 
pative effect is accounted by conventional bulk viscous effect. 
In FRW universe the bulk viscosity can be modelled within 
the framework of non-equilibrium thermodynamics proposed 
by Israel and Stewart [31]. According to this theory the bulk 
viscous pressure tt follows the transport equation 



-3nH 



[3H 



T 
T 



T 
T 



(21) 



where the positive definite quantity r] stands for the coefficient 
of bulk viscosity, T is the temperature of the fluid and r is 
the relaxation time associated with the dissipative effect i.e, 
the time taken by the system to reach equilibrium state if the 
dissipative effect is suddenly switched off. Considering the 
divergence term in the square bracket to be small i.e, (A^r) 
to be constant, the equation can be approximated to a simpler 
form 



tix = —3nH 



(22) 



In literature this is commonly described as a truncated version 
of the full nonequilibrium thermodynamics. The viscous ef- 
fects are assumed to be not so large as observation seems to 
rule out huge entropy production on large scales p2]]. Usu- 
ally t is expressed as 2 so as to ensure that the viscous signal 
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does not exceed the speed of light \ p3\ j and also (tH) = v, 
where v > 1 for a consistent hydrodynamical description of 
the fluid fl34|]. With these two assumption equation ( [22"| ) be- 
comes 



3/5^ 



(23) 



In a very recent work Chimento et al [ ]26| ] have shown that a 
mixture of minimally coupled self interacting scalar field and 
a perfect fluid is unable to drive accelerated expansion and 
solve the cosmic coincidence problem at the same time, while 
the mixture of a dissipative CDM with bulk viscosity along 
with minimally coupled self interacting scalar field can suc- 
cessfully drive the accelerated expansion and solve the cosmic 
coincidence problem simultaneously. 

An effective negative pressure can also be created from cos- 
mic antifriction which is closely related to particle production 
out of gravitational field. In a recent paper Zimdahl et al [ |28| ] 
have shown that one can have a negative it if there exists a 
particle number nonconserving interaction inside matter. This 
may happen due to particle production out of gravitational 
field. In this case, the matter is of course not a dissipative 
fluid, but a perfect fluid with varying particle number. Though 
substantial particle production is an event that occurs in early 
universe, Zimdahl et al have shown that extremely small par- 
ticle production rate can also cause sufficiently negative 7r to 
violate strong energy condition. 

We do not apriori assume any specific model for this dis- 
sipative effect, rather only assume the existence of a negative 
7T. For a similar kind of evolution of the scale factor and the 
scalar field given by equation (|l2|), the energy density for the 
fluid, with negative pressure is also given by equations ( pj[ ) 
and (jl4|). From equations (|o|) and (jl3|), one can easily find 
that 



(2-P- 3a)c/> 
t 13 



2a + (3(l + a)-(3 2 (l+Lu) 



2-(3 



.0-2 



(24) 



From equation ( [24] ) one can easily check that to have a nega- 
tive 7r, one should have 3a > 2 — (3 which essentially means 
a > §. This suits the condition for acceleration as an a is 
needed to be greater than 1 for that. One can also check that 
for this set of solutions given by @, C|l3j> and @ and from 
equations (|23"|), the condition v > 1 holds provided 2 — (3 > 0, 
which is very much true in our case. This is important for the 
hydrodynamical description if the CDM is assumed to be a 
conventional viscous fluid. 

To have a clear picture of the expansion of the universe and 
the missing energy, we further study the energy density and 
pressure of the geometric scalar field. The expressions for 
the energy density and the pressure of the scalar field can be 
derived from the field equations (J^) and ([))) to be 



w </> 2 V R 



(25) 



and 



2 ~~6 



V 
~2 



4* 



In case of power law expansion (|12j) and potential like i 
energy density of the BD field becomes 



Ptj> = <i 3c i 

and pressure of the BD field is 

a (n 5 
Pv = ~2 V 2 



- a 



2 

—uj - 
3 



(26) 



i the 



(27) 



(28) 



The positivity condition for the scalar energy density demands 



w >-(3a+-) 



(29) 



which eventually restricts u> beyond some lower value. So 
essentially the two positivity energy conditions, (|l8]) and (j2^), 
limit the range of uj within —(3a + 1) < u> < 3a— |. Clearly 
a barotropic relation (p^ = is followed by the scalar 

field, where the adiabatic index 7^ is given by 



14> = - 1 



(a + 



f (3c 



(30) 



The range of 7^ that agrees with the observational datas and 
describes the current acceleration for the universe well, is 
—0.6 > 70 > —1. One can adjust the value of a and ui so 
as to get the required value of 7^. We now recast equation (^|) 
in the form 



f2 m + fi^ — 1 



(31) 



where the density parameters for matter fi m and scalar field 
f2i* are defined to be (see ref [B5]]) 



and Cls = 



The expression for density parameters at present epoch are 



(32) 



^mO — 



Po = 1 _ J_/ 5x 
3Hg4> 2 6a [UJ 2 } 



and 



n 



(p0 



P<t>0 



11, 5. 
i, -0 + 7^+2' 



3.ff^o 2 6a 



(33) 



(34) 



Like 70, the value of 57 m0 that suits best the luminosity 
distance-redshift data for type la supernovae is Q m o — 0.28 
and in a similar fashion like 7^, one can adjust a and u> value 
to get the required value of Q m o that tally with the observa- 
tion. 
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So far a is restricted only by the deceleration parameter 
that it should be greater than 1 for the universe to accelerate. 
And the positivity energy conditions limits oj within the range 
— (3a + |) < o-i < 3a — |. So if a small value is chosen for 
a (~ 1), oj is also restricted accordingly and a suitable choice 
of both parameters can be made to find the allowable range 
of 70 and f2 m o that matches observation. Kaplinghat et al 
[ |36| ] and others [ f37| ] have pointed out that for power law cos- 
mologies, high redshift data and present age of the universe 
restricts a to a value « 1. In a very recent investigation we 
have found that the best fit value of a with the SN la data for 
power law cosmology is approximately 1.25 J38| ]. Hence this 
small value for a restricts u> also to be small. But this squarely 
contradicts the solar system bound on oj(> 600). To accom- 
modate such large value of oj, a should be large. But for the 
large values of a the universe accelerates faster ( almost like 
de-sitter expansion) and 70 asymptotically approaches —1. It 
is quite unlikely that the universe presently accelerates in such 
a high power law fashion and such a values do not match with 
present observation. We would discuss more about this point 
in the discussion section. 

We wish to find the range of the parameters a and oj al- 
lowable in our model , that suits the permissible range of 
70 and tt m o of the quintessence proposals. In figure 1, we 
have shown the allowed region in the (a, oj) parameter space 
(shaded portion in the figure) for the specified range of 70 
(—0.6 > 70 > —0.8) and O m0 (0.5 > fi m o > 0.3), where we 
have assumed a is small. It can be noticed that the allowable 
range of both the parameters a and to obeys the constraints, 
(p^)and (p9|), imposed on them by physical conditions. 
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FIG. 1. ujvsa for -0.6 > 70 > -0.8 and 0.5 > fi m0 > 0.3 

A point to note here is that in BD theory the gravitational 
coupling G varies inversely with the scalar field </>. At present 
time <fi approaches a constant value </>o, the inverse of which 
gives the present Newtonian constant Gn- In the weak field 
limit the present Newtonian coupling and the asymptotic value 
of <f> is related by 



2uj + A 1 

2u> + 3 <fio 



(35) 



The present day variation of the gravitational coupling G, is 

§1° = Ja^n, where Ho (= j^) is the Hubble parameter at 
present. For any value of a that allows acceleration this rate 
is < 10~ 10 per year 

Another important point to mention here is that time vari- 
ation of G does not directly affect the nuclear process of the 
early universe. But the expansion rate of the universe in this 
type of theory do influence the primeval nucleosynthesis [p0|]. 
A fixed value for the parameter a > 1 for all epochs implies 
that universe is always accelerating which seriously contra- 
dicts the nucleosynthesis scenario. One way to avoid such 
problem is to consider oj as a function of the scalar field <f>. In 
a recent work Banerjee and Pavon [[HJ have shown that with 
oj(4>) one can have a decelerating radiation dominated era in 
the early time and accelerated matter dominated era in the late 
time. But in their case also, oj has to be small asymptotically 
to have a late time acceleration for the universe. 

To analyse the nature of acceleration and our ansatz more 
critically, it is interesting to match different cosmological pa- 
rameters with observations. We intend to find the age of the 
universe and the luminosity distance-redshift relation compat- 
ible with our model, probing the background dynamics, that 
could differentiate between different types of universe. 

Since one of the main incentive for reconsidering the in- 
troduction of cosmological constant was the age of the uni- 
verse, we first consider the age of the universe suggested in 
our model and the constraints imposed on it by observations. 
Equation (Jfy can also be presented as 
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FIG. 2. Age (t) vs redshift z for different a values 
where z is the redshift defined by 
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From equation ( pq ) we find the age of the observable uni- 
verse for a given redshift z, is 



t -t 



Ho(^mO + ^00 ) 2 



1 



1 



(1 



(38) 



Of course for t — 0, i.e, for infinite redshift, the age of the 
universe is to = 



H 



An old object observed at a certain redshift selects all mod- 
els with at least that age at that given redshift. In that respect, 
several age constraints have recently appeared in the literature 
[p^]. For example, the age of the radio galaxy 53W091 ob- 
served at a redshift z — 1.55 puts a lower bound of 3.5 Gyrs 
at that redshift. The quasar observed at z = 3.62 sets a lower 
bound of 1.3 Gyrs. In figure 2 we present a plot of the age of 
the universe as a function of redshift for various values of a. 
Taking into account the range of a, prescribed by figure ([!]), 
our universe has an age limit of to > 15.5Gyrs. From the 
figure it can be seen that all the universes with that minimum 
age limit , except that of Einstein-desitter one, always accom- 
modate these constraints. Recently Pont et al [[l3| estimated 
the age of the universe to be 14 ± 2 Gyrs, which is in excellent 
agreement with our result. 

Now we would like to trace the change of luminosity dis- 
tance with respect to the redshift in our model so as to com- 
pare it with the present data available. The result that reveal 
the so called acceleration of the universe was the observa- 
tion of the luminosity distance as a function of redshift for 
type la supernovae, which is believed to be a standard candle. 
From almost 60 redshifts, 42 high redshift data obtained by 
Supernova Cosmology Project and 18 low redshift observed 
by Calan Tololo Supernova Survey, favour a universe with 
positive cosmological constant. Assuming flatness in con- 
text of general relativity, the best fit for these data occurs for 
Sl m o = 0.28 and f^o = 0.72. Optical astronomers measure 
luminosities in logarithmic units, called magnitudes, given by 



m(z) = M + 5 log d L + 25 



(39) 



where M. is the absolute magnitude and d^ is luminosity dis- 
tance defined by 



d L = R(t ){i + z) n 

for an event at r = n at time t = t\. 
According to our ansatz the expression for d^ is 



d L (z) 



where F(z) = 



(1 + *) 



Ho(fl m o + fl<f>o) 2 Jo 



F(zf)dzf 



(40) 



(41) 



(l+z)- 



In figure 3 we have plotted this luminosity distance versus 
redshift for different values of a. We see that for different 
a values, the di is practically same for lower redshifts upto 
z ~ 0.4. At redshifts z > 0.4 the curves are separated, but 
is not distinctly separate to discriminate and rule out different 



types of the models. Therefore, high accuracy measurements 
with uncertainties at percentage level are needed in order to 
cleanly distinguish the models and the need to go to redshifts 
sensibly higher than 1 is evident. In this respect it is very 
much relevant to mention that Supernova Acceleration Probe 
(SNAP) is planned to make measurements with an accuracy 
at percentage level upto redshifts z ~ 1.7 
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III. DISCUSSION 

This work investigates the possibility of getting an accel- 
erated universe in context of BD theory with a 4 potential 
having a time dependent mass squared term and a matter field. 
In this work we have not used quintessence field to trace the 
missing energy. The BD scalar field , which is a geometric 
scalar field, plays the role of dynamical A and provides that 
missing energy. It is found that for a simplistic approach of 
power law expansion (~ t a ) a perfect fluid kind of matter 
(both pressureless and with pressure) cannot support a late 
time acceleration of the universe, if the scalar field have po- 
tential given by equation (^). But a matter with a dissipative 
effect can provide the acceleration that agrees with the ob- 
servational data sets. The dissipative effect accounted by the 
negative pressure can be modelled in two ways according to 
the recent investigations J2^^6| ]. Particle production out of 
gravitational field can give rise to negative pressure while en- 
ergy can also be dissipated by bulk viscous effect between the 
CDM particles. In this work though no particular model is 
considered for the origin of the negative pressure, it is found 
that v > 1. This is important for hydrodynamical descrip- 
tion if the CDM is assumed to be a conventional dissipative 
fluid. We have also calculated different parameters like the 
time variation of gravitational coupling, age of the universe 
and the luminosity-distance redshift relation. All of these cos- 
mological parameters agrees quite well with the recent obser- 
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vations. 

The accelerated solution depends crucially upon two pa- 
rameters: a and BD parameter u, both of which are con- 
strained by different physical conditions. Different combi- 
nations of a and to can produce the required values for 7^ 
(—0.6 > 70 > —1) and VL m (~ 0.3) that tallies with present 
observation suggesting acceleration. Small a values restricts 
u> to small negative values and support the late time accelera- 
tion scenario quite successfully. The cosmological parameters 
calculated with the small value of a agrees with observations 
quite well. Many references [ p6[[35| , [37| ] are available in the lit- 
erature where it has been shown that this value should be very 
close to 1 to be consistent with observation. In one of our re- 
cent work [ f38t ] it is shown that for simple power law expansion 
(~ t a ) of the universe, the best fit value of a with SNIa data is 
approximately 1.25, and due to jl^ ) and (^9|), this will restrict 
uj to small value. But this clearly contradicts the solar system 
limit uj > 600. One should note that in our model, large value 
of id consistent with solar system limit, is not restricted either 
by physical conditions such as positivity of energy density or 
by the requirements of specific ranges for 7^ or Q m consistent 
with the observations. Only the fact, that a is not large, is the 
prediction of the datas obtained so far and this constrains id to 
small value. If future observations predicts large value for a 
then that can also be accommodated in our model with a large 
value for u>. 



lem is to consider u> to be function of the scalar field cf> [Sjj]. A 
choice of id (polynomial function of <p) can give a decelerating 
radiation era as well as accelerating matter dominated era. But 
then also id asymptotically acquires a small negative value for 
an accelerating universe at late time. In most of the investiga- 
tions done in scalar tensor theory [ ^4] , ^8|p4[ | such conclusions 
have been arrived. Only exception, so far in our knowledge, 
is that done by Bertolami and Martins [22], where in Brans 
Dicke cosmology with a <fi 2 potential, the solution of an accel- 
erated universe is obtained with large |a>|. But there the pos- 
itive energy condition on both matter and the the scalar field 
have not been considered. There are also other evidences in 
literature where small |o;| has been supported. In the extended 
inflationary model by La and Steinhardt [Q], the required 
value for id is 20. The structure formation in scalar tensor the- 



ory also contradicts the solar system bound on id [ 45 1. Thus 
the problem seems to appear in different scales (astronomical 
and cosmological). The theory has been tested by experiments 
so far only in the astronomical scales and no experiment had 
been done in cosmological scale as yet. And so the problem 
occurs in finding the compatibility between astronomical ob- 
servation and cosmological requirements. Considering id to be 
a variable to have both decelerating and accelerating phases at 
different epochs, while large id values occur due to local inho- 
mogeneities in astronomical scale to satisfy the solar system 
bound, can be a complete investigation and may give a satis- 
factory answer to this question. 



One should also note that, for a < 1.33 which is consistent 
with the present data, the dissipative pressure is not sufficient 
to drive the acceleration alone, and the BD scalar field along 
with the dissipative pressure in the CDM drives the accelera- 
tion whereas the BD scalar field plays the role of the missing 
components of the universe. 

It is also expected that after the phase transition, one ex- 
pects 4> to roll down from = taking different values in dif- 
ferent directions causing large scale inhomogeneities. But it 
was argued by McDonald [|3(j] that after this recent phase tran- 
sition the universe will be filled by non topological objects like 
"axiom" whose radius is much smaller than lOAlpc caus- 
ing these axions to behave like a smooth dark energy compo- 
nents. 

The range of id, we obtained in our calculations in order 
to fulfil different physical conditions and also to have correct 
range of values for 7^ , a, and f^mjis consistent with the range 
obtained by other authors [ ^J^ffl] ] . Also, it can be seen in 
the figure 3 that for different values of id and a, it is difficult 
to distinguish between models upto z ~ 1. Hence our model 
is not very much fine tuned as far as the parameter a and id 
are concerned. 

It is also important to note that as a remains constant 
throughout the age of the universe. This essentially means 
that the universe always accelerates for a > 1 , which se- 
riously contradicts the primeval nucleosynthesis scenario. As 
we have mentioned earlier that one way to overcome this prob- 
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